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We estimate the ratios of bare astrophysical S-factors at zero incident energy for proton and 
deuteron induced reactions in a model which assumes a compound nucleus formation probability 
plus a statistical decay. The obtained ratios agree well with available experimental values, as far 
as the reactions which have dominant s-wave entrance channel components are investigated. Due 
to its simplicity the model could be used as a guidance for predictions on reactions which have not 
been investigated yet. 



PACS numbers: 

The nuclear fusion cross sections of proton and 
deuteron induced reactions at low energies are of par- 
ticular interest from the points of view of the stellar nu- 
cleosynthesis and the nuclear energy production. These 
cross sections are measured at laboratory energies and 
extrapolated to thermal energies 0, 0] , because of their 
small values at such low energies. This extrapolation is 
done by introducing the astrophysical S-factor: 

S(E) = a(E)Ee 2 ^ E \ (1) 

where <r(E) is the reaction cross section at the incident 

center-of-mass energy E and r)(E) = Z^Zpa-J Zf, 
Zp, [i denoting the atomic numbers and the reduced 
mass of the target and the projectile, a and c are the 
fine-structure constant and the speed of light, respec- 
tively. The exponential term in the equation represents 
the Coulomb barrier penetrability. Since one has fac- 
tored out the strong energy dependence of <j(E) due to 
the barrier penetrability, the S-factor could be approx- 
imated by a smooth energy dependence in the absence 
of low-energy resonance. A lot of effort has been made 
to extract the S-factor in the low energy region inpar- 
ticular for transfer reactions experimentally [1, 0, [a [|| . 
Although to know the low energy S-factor for radiative 
capture reactions, still one needs the extrapolation with 
theoretical formulae 0). Owing to this experimental ef- 
fort one can, in principle, determine the S-factor in the 
energy region of astrophysical interest directly. However 
at such a low energy it is known that electrons around the 
target nucleus have an effect on the fusion cross section. 
In contrast, in the stellar nucleosynthesis, nuclei are al- 
most fully ionized and are surrounded by the plasma elec- 
trons. The nuclear reactions in such a circumstance are 
affected by a different mechanism of the plasma electron 
screening [8, 9]. Hence the screening effects of the bound 
electrons should be removed from the S-factor data, in or- 
der to asses the reaction rate in the stellar site correctly. 
The enhancement by the bound electrons is discussed 
in terms of a constant potential shift (screening potential 
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In this paper we would like to investigate in detail the 
physical origins of the bare S-factor. We propose visual- 
izing the fusion mechanism under the assumption of the 
compound nucleus (CN) formation as a two step process: 

STEP 1) The nuclear fusion under the Coulomb barrier 
occurs with a given probability which depends crucially 
on nuclear effects (the nuclear potential, nuclear sizes, 
an excitation of collective modes etc.). A CN is formed 
at a given excitation energy. Since we are interested in 
determining the S-factor in the limit of zero bombarding 
energies, the excitation energy is essentially given by the 
Q-value of the reaction which leads to the CN. 

STEP 2) The CN decays into a given channel. We 
assume this decay process to be statistical and make use 
of the Weisskopf model [2 EH to evaluate its probability. 

Low energy behaviors of the S-factors have been stud- 
ied by means of the R- matrix theory [l6|, U3, EH , which 
is also based on a description of the nuclear reaction 
through the intermediate CN. In the nuclear reactions in 
general, it is known that this CN process coexists with 
the direct process mechanism [l9j . Provided that the 
energy regions of our interest in this paper is low, we 
presume that the / = partial wave component would 
be dominant for most of the reactions investigated. This 
fact suggests that most of low energy reactions are pos- 
sibly to be described by using only the CN mechanism. 
Certainly there are some reactions which do not have 
I =0 component, these cases will be discussed separately. 
The situation of our interest is displayed in Fig. [TJ where 
the colliding ions at a given beam energy E in the center 
of mass (CM) system must overcome the Coulomb bar- 
rier in order to fuse. In general, the cross section for the 
sub-barrier fusion at energy E and angular momentum I 
is given by [20( | : 

a(E) = TTh 2 /(2 t iE) 5^(21 + l)Ti(E), (2) 

where Ti(E) is the tunneling probability. We take the 
limit of E — y and keep only I = term of the sum in 
Eq. |(5J) in the following discussion. However one can eas- 
ily extend the case where the centrifugal potential term 
is included [2l| and later we will only mention the mod- 
ification by the I > 1 terms to the I = case. We point 
out that Ti(E) is not the pure Coulomb penetrability but 
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FIG. 1: A schematic feature of the interaction between col- 
liding nuclei. 



the penetrability of the Coulomb + centrifugal + nuclear 
potentials. The tunneling probability Tq(E) can be writ- 
ten in terms of the action A in the low incident energy 
limit USUI: 



T (E) ~exp 









, A= f 




Jr 



.(0) 



dry/2p[V c {r) - E] 



(3) 

where V c (r), Rn and are the Coulomb potential, the 
inner and the outer classical turning points, respectively, 
as they are shown in Fig. [TJ In particular, if we consider 
the pure Coulomb penetrability, i.e. we take the limit of 

Rn — ► 0, A becomes Ac = ZpZpe 2 \J~^§-: in the Gamow 
limit. As a consequence Eq. ^ reduces to the Gamow 
factor which is exactly the inverse of the exponential term 
of Eq. (QJ. Comparing eqs. ^ and ^ we derive an S- 
factor in the Gamow limit: Sq — 7rft 2 /(2/i). It means 
that the S factor would be independent of the energy and 
of the exit channels, if the process would be dominated 
by Coulomb. However the actual value of the S-factor 
depends on the colliding ions and on the exit channels. 
For instance let us consider two reactions 10 B(<i, a) 8 Be 
(S exp = 102.7 MeV b) [HJ and n B(p, a) 8 Be (S exp = 
208.0 MeV b) [23[ which have identical exit channels but 
different entrance channels. Our simple formula gives 
Sq = 0.39 and 0.71 MeVb, respectively, which disagree 
with both experiments. This simple comparison suggests 
that the nuclear effects are really dominant and should 
be taken into account. A modification due to the nuclear 
effects can be easily included assuming that the nuclear 
potential is given by a square well of radius Rn indicated 
by the dotted line in Fig. [T] [2l[ . Here we will assume the 
empirical formula: Rn — a n RNo, where Rno = 1.2 x 



(4 /3 - 



1/3 

-Arp ) [fm] is the sum of the radii of the projectile 



and the target, denoting the mass numbers of the target 
and the projectile nuclei by At and Ap, respectively. a n 
is a parameter which takes into account the fact that 



nuclei have a diffuse potential and during the reaction 
some collective degrees of freedom might be excited which 
could modify the effective radii of the ions. 

In the approximation considered, the action is written 
down in terms of Ac as: 



A = Ar 



dr^2p [V c (r)-E], 



(4) 



where the integral of the right hand side is the action of 
the pure Coulomb potential from the origin to the nuclear 
interaction radius. This action integral is approximated 
by 2-\/2/iZyZpe 2 i?jv, in the limit of E — > and therefore 
the S-factor modified by nuclear effects is written by, 



(5) 



A simple inspection of Eq. ([5]) using standard values is 
not able to reproduce the experimental S-factors for the 
two reactions discussed above. This implies that the 
modification is not enough, in fact we have only included 
STEP 1) of our proposed scenario. We need now to mul- 
tiply Sno by the probability of getting a certain exit 
channel from the CN U(A* -»• c + d). The Weisskopf 
modelfllj gives an expression for the probability to ob- 
tain a certain exit channel i as a function of the cross 
section of the reverse process af bs : 



n i =g i (Tl+2m i T ki ) 



CN 



P(E* CN ) 



(6) 



One can, therefore, express the probability n as 
ELj/5^nj, where the sum in the denominator is taken 
over all possible exit channels. In Eq. (O <fa denotes the 
number of states for the spin of the particle considered. 
(<7; = 2sj + 1, where Sj is the spin for proton, neutron 
and alpha particle, gi — 2 for 7 [la].) and m, are the 
kinetic energy and the mass, respectively, of the light- 
est reaction product in the exit channel. pi(ED N — 
and p(Eq N ) are the level densities of the heavier reaction 
product and of the CN, respectively, where E^ N is the 
excitation energy of the CN. Introducing the entropy of 
the nuclei, the ratio of pi(Ec N — T^J to p(E c \ N ) can be 
written as 



Pi( E CN ~ T kj) 

p(E* C n) 



ex p( — y>' 



(7) 



where T is the temperature of the CN, in the u nit of cn- 
ergy. The CN temperature is given by T — yQcN / a , 
Qcn is the Q- value for the CN formation and a ~ 
Acn/&-0 [MeV -1 ] is the level density parameter [3]. 
Here we have denoted the mass number of the CN to 
be Aqn and have assumed that all the excitation energy 
is given by the Q-values (i.e.E — > 0). This is the re- 
quired STEP 2) of our model. We can, therefore, define 
a S-factor of the I — entrance channel as: 



s G e* y/2 ' J - ZTZpRN n. 



(8) 
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The factor depends on the exit channel as expected. In 
passing we mention that the S'-factor of the reactions 
which have I > 1 partial wave components in addition to 
the I = component can be written by [2l[ 



S N = s Ge tV2^Ffl7 



i + ym + i)e 



i>i 



n. 



(9) 

One can, in principle, obtain the absolute value of the 
S-factor from Eq. (JSj) , however to determine the level 
density p(E CN ) one needs some numerical calculations. 
We will, therefore, discuss limiting cases where the model 
gives an analytical solution, in order to get a feeling of the 
proposed approach. The first case is given by the ratios 
of S-factors (in the limit of zero beam energy) where the 
CN formed is the same but two different reaction chan- 
nels are studied (different isospins in the exit channel), 
while the second is the opposite i.e. different isospin in 
the entrance channel and the same decay mode. In this 
way, by taking ratios, screening effects should exactly 
cancel in both cases. 

When the two reactions with an identical entrance 
channel are considered, its S-factor ratio R is given by: 



R = 



9 An, 



2mjT kj 



,-AQ/T 



g i (Tl+2m i T ki )a. 



abs 



(10) 



where refers to the decay channels, and Tfe., 4) = 

A+A Q?'(0' reaction Q- value for each channel and 
its difference is denoted by QjU) and AQ = Qj — Q il 
respectively. We write the mass numbers of particles in 
the exit channels as A\ for the lightest particle and 
for the rest. Thus in the defined ratio the 'entrance chan- 
nel' effects drop out. In the case of the radiative capture 
the equation above can be easily modified using the rel- 
ativistic Weisskopf equation given above, eq.©. The 
absorption cross sections are taken to be geometrical but 
modified by a Coulomb barrier: 



abs 



TzR 2 N {l-V c (R N )IQj), 



(11) 



where Rn is the sum of the radii of the nuclei in the 
exit channel j. For the photo- absorption cross section 
we make use of the parameterization given in [24j |. 



abs 



75A 



CN- 



(T fc7 - 2.22 6) 
T 3 



3/2 



[mb], 



(12) 



where T kl = Qcn- This parameterization holds well for 
relatively high 7 energies (T^ 7 > 23MeV). We will make 
use of it, if there is no experimental data available at 
lower energies. 

Within this simple approach we have estimated the 
S-factor ratios for the proton and deuteron induced re- 
actions with positive Q-values on the nuclei D, 6 Li, 7 Li, 
9 Be, 10 B and 11 B. All the investigated reactions are listed 
in Table U and HT] We take the ratios of the 2nd to the 
3rd column reactions. Their S-factor ratios determined 



TABLE I: Ratios of S-factors for the transfer reactions with 
identical entrance channels and different exit channels. 



Reactionl 
So(MeVb) 


Acn 


Reaction2 
So(MeVb) 


Ratio 
(data) 


Ratio 
(calc.) 


d(d, p)t 
0.0571 [17J 


4 


d(d,n) 3 He 
0.0524 [17] 


1.09 


0.97 


6 Li(d,n) 7 Be 
25.3 


8 


6 Li(d,p) 7 Li 
26.6 [25] 


0.95 [25] 


1.26 


6 Li(d, a) 4 He 
23.1 [26] 




6 Li(d,p) 7 Li 
26.6 [25] 


0.87 


0.37 


9 Be(p, d) 8 Be 
18. [1] 


10 


9 Be(p,a) 6 Li 
18. [Tj 


1.0 


0.95 


9 Be(d,p) 10 Be 
15.2 [2, 22J 

9 Be(d, t ) 8 Be 
5.58 [2, 22] 


11 


9 Be(d,a) 7 Li 
20.8 [2, 22] 
9 Be(d,a) 7 Li 
20.8 [2, 22] 


0.73 

0.26° 


0.79 
2.3 


10 B(d,a) 8 Be 
102.79 [2, 22] 


12 


156.1 [2, 22] 


0.66 


0.42 


n B(d,a) 9 Be 
93.2 [2, 22J 


13 


n B(d,p) 12 B 
105. [2, 22] 


0.89 


1.72 



a Experimental data for the reaction 9 Be(ei, io) 8 Be is 
measured only for the 8 Be ground state. 



from experimental data are shown in the 4th column and 
are compared with the ratios from Eq. (|10p in the last 
column. In Table HT1 we made a list of the ratios of the ra- 
diative capture reactions to the transfer reactions. From 
the tables we see a general good agreement of our esti- 
mations to data. In the Acn = 11 case in the Table U 
the model gives a ratio much larger than data for the 
reactions 9 Be(d, t) 8 Be/ 9 Be(d,p) 10 Be, however we note 
that the experimental value for the tritium case refers to 
the to channel case only. According to our approach, we 
would estimate an astrophysical factor for the t-channel 
at least a factor 10 than what is quoted in [22j]. For the 
same CN case but for the 7 channel in the Table HT1 we 
reproduce the ratio within a factor of two. This is sur- 
prising because this reaction is dominated by a resonance 
at about 9 KeV and we have estimated the experimen- 
tal ratios from the extrapolation of ref . [l| . Evidently the 
effects of the resonance cancel out by taking ratios. The 
only discordance between a calculated ratio and data is 
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TABLE II: Ratios of S-factors for the radiative capture reac- 
tions to the transfer reactions with identical entrance chan- 
nels. 



Acn 








Reactionl 


Reaction2 


Ratio 


Ratio 


S'o(MeVb) 


So(MeVb) 


(data) 


(calc.) 



d(d,i)t <f(d,n) 3 He 
0.005 x 10" 6 [1] 0.0524 [17] 9.5xl0" 8 7.9xl0 -8 a 



7 



3 Li(p,7) 7 Be 6 Li(p,«) 3 He 



1.07xl0 -4 [1] 



3. [1] 



3.6x10" 



9.6x10" 



8 



7 Li(p, 7 a) 4 He 7 Li(p,a) 4 He 



0.0015 [1] 



0.074 [1] 



10 



9 Be(p,7) 10 B 9 Be(p,d) 8 Be 
0.001 [1] 18. [1] 

9 Be(p, 7 ) 10 B 9 Be(p,a) 6 Li 
0.001 [1] 18. [1] 
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10 B( P ,7) U C 
0.00170 [1] 



10 B(p,a) 7 Be 
70. [1] 



2.0xl0" 2 1.4X10" 36 



5.6xl0 -5 6.2xl0 -5 
5.6xl0 -5 5.9xl0~ 5 

2.4xl0" 5 5.2xl0" 5 



12 



l B(p,7) 12 C 11 B(p,a) 8 Be 



0.0035 [1] 



200.0 [1] 1.8 xlO" 5 2.4x10" 



16 

N(p, 7 ) 16 15 N(p,a) 12 C 



0.064 [1] 



65. [1] 



9.8x10" 



1.8x10" 



a He(7,d)D has the photo disintegration cross section 3.2 [fi 
b] at the peak of T fc7 =29MeV [27| 
b In the reaction 7 Li(p,7a) 4 He we have assumed that 7 
kinetic energy is equivalent to the CN temperature. 



seen in the case of Aqn = 8 in Table HT1 This disagree- 
ment of the obtained ratio arises from the fact that the 
reaction Li(p, ev) 4 He does not have I —0 partial wave in 
the entrance channel, because of the parity conservation 
in the exit channel. This reaction is restricted to have 
odd number partial waves in the entrance channel. As 
we mentioned in the introductory part, our model de- 
veloped in this paper assumes s-wave entrance channel 
for simplicity. Therefore our model fails to estimate the 
ratio fot this case properly. We stress that since the en- 
trance channel is the same when we take the ratios, an 
eventual screening should cancel out exactly. Thus our 
method should give some constraint on extracted bare 
astrophysical factors. 



TABLE III: Ratios of S-factors for transfer reactions with 
similar Q- values. 





Acn 






Reactionl 


Reaction2 


Ratio 


Ratio 


So(MeVb) 


So(MeVb) 


(data) 


(calc.) 


D(£,n) 4 He 


5 

3 He(d,p) 4 He 






10. [1] 


6.3 [18] 


1.6 


0.42 



Another interesting example is 

D(i, n) 4 He/D( 3 He, p) 4 He: the ratios of reactions 
with similar Q- values. Using eqs. §5Mj§ and assuming 
that the correction coming from the CN decay is given 
by their respective Q-values and temperatures, the ratio 
of the S-factor is expressed as [2lj |: 



R 



fJ-j exp (j- v / 2fi l Z Tz Z Pi R 



IX 



(13) 



where i and j refer to the different entrance channels. 
For the present case, both reactions have the only main 
exit channel, so that the statistical terms are simplified. 
The ratio obtained by Eq. (fTB")) , which is shown in the 5th 
column in the Table IIII1 is smaller than data by a factor 
of 3. If we substitute a n . by 2.5, we can reproduce the 
experimental ratio. 

Last examples we consider the ratios of reactions with 
different isospins in the entrance channel but the same 
exit channel. We consider two ratios for the reactions 
listed in Table IIVI All the possible exit channels (i.e., 
n, p, d, i, 3 He and a) are included to evaluate Eq. (fT3"|) . In 
the 5th column the ratios obtained by Eq. (|13p are shown. 
For the first ratio Eq. (fl~3|) underestimates the data by a 
factor of 1000, while for the second ratio Eq. (fT3")) gives a 
value in agreement with the data. The disagreement of 
the obtained ratio for the reactions which include lithium 
isotopes again is caused by the absence of the I =0 partial 
wave in the entrance channel of the reaction 7 Li(p, a) 4 He. 
By contrast, in the case of boron involved reactions both 
reactions have the I =0 component 18, 28] and as a mat- 
ter of fact our model gives reasonable ratio for these two 
reactions. If we use the the I — 1 term of the formula © 
for the reaction 7 Li(p, a) 4 He, we obtain 9.0, instead of 
0.25, as the ratio. We can reproduce the ratio by substi- 
tuting a n . by 3.3 in the former case. It means that the 
nuclear interaction radius is larger for the entrance chan- 



18], 



nel 6 Li+d, which has a dominant s-wave component 
of a factor of 3.3 with respect to the reaction 7 Li+p. 

These two last cases suggest that nuclear effects play 
an important role and should be considered in detail. 

In conclusion, we have discussed the possible physi- 
cal origins of the astrophysical S-factors. We have shown 
that the essential contributions to the S-factor come from 
the Coulomb and nuclear penetrability of the colliding 
nuclei to form a CN and from its decay in a given chan- 
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TABLE IV: Ratios of S-factors for transfer reactions with 
different entrance channels and identical exit channels. 



Acn 

Reaction 1 
So(MeVb) 


Reaction2 
S'o(MeVb) 


Ratio 
(data) 


Ratio 
(calc.) 


8 

6 Li(d, a) 4 He 


7 Li(p, a) 4 He 






16.9 [4] 


0.055 [5J 


307. 


0.25 


12 








10 B(d,a) 8 Be 


n B(p,a) 8 Be 






102.7 [2, 22] 


208.0 [23] 


0.49 


0.46 



nel. We have assumed that this decay process is described 
statistically. By taking suitable ratios of astrophysical 
factors, one can treat those two ingredients separately 
and compare the model predictions easily to the exper- 
imental data. In the case where the 'entrance channels' 
effects are cancelled out in the ratio, we have shown that 



our model can reproduce the experimental data without 
any fitting parameter. When we consider the cases where 
the same CN is formed in different reactions, the nuclear 
effects become important. When we take the S-factor ra- 
tio of a radiative capture reaction to a transfer reaction, 
this approach can predict the wide reduction of the yield 
in the radiative capture. Our model demonstrates clearly 
that this reduction arises from the ratio of the 7 kinetic 
energy and the ion mass. Our estimations have been per- 
formed in the limit of low incident energies, so that our 
results refer to ratios of bare astrophysical factors at zero 
incident energy. At such low energies experimental data 
of S-factors are known to be influenced by surrounding 
electrons. The procedure of taking the ratios of two reac- 
tions with an identical entrance channel allows us to give 
some constraint on extracted bare astrophysical factors. 
Having acquired some confidence in our model, we plan 
to make use of it to reproduce the experimental data of 
S-factors as a function of the beam energy. 
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